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The Bogoliubov method for the excitation spectrum of 
a Bose-condensed gas is generalized to apply to a gas with 
an exact large number N of particles. This generalization 
yields a description of the Schrodinger picture field opera- 
tors as the product of an annihilation operator A for the to- 
tal number of particles and the sum of a "condensate wave- 
function" £(x) and a phonon field operator x( x ) m the form 
i/)(x) w j4{^(x) + x( x )/V^V} when the field operator acts 
on the N particle subspace. It is then possible to expand 
the Hamiltonian in decreasing powers of v~N, an thus obtain 
solutions for eigenvalues and eigenstates as an asymptotic ex- 
pansion of the same kind. It is also possible to compute all 
matrix elements of field operators between states of different 
N. 

The excitation spectrum can be obtained by essentially the 
same method as Bogoliubov only if £(x) is a solution of the 
time independent Gross-Pitaevskii equation for N particles 
and any chemical potential fi which yields a valid and sta- 
ble solution of the Gross-Pitaevskii equation. The treatment 
within a subspace of fixed N is identical in form to that usu- 
ally used, but the interpretation of the operators is slightly 
different. 

A time-dependent generalization is then made, yielding an 
asymptotic expansion in decreasing powers of \J~N for the 
equations of motion. In this expansion the condensate wave- 
function has the time-dependent form £(x, £), and the con- 
dition for the validity of the expansion is that f(x, £) satis- 
fies the time-dependent Gross-Pitaevskii equation d^/dt = 
-(h 2 /2m)\/ 2 ^ + V^ + Nu\^\ 2 C 

The physics is then described in a kind of interaction pic- 
ture, called the condensate picture, in which the phonon op- 
erator can be expressed as x( x >*) = 5Zfc £ fe ( x > t) a k, where 
the operators are time independent annihilation opera- 
tors, and the state-vector has a time evolution described by 
a Schrodinger equation in which the Hamiltonian is a time- 
dependent quadratic form in the phonon creation and an- 
nihilation operators, whose coefficients are explicitly deter- 
mined in terms of the time-dependent condensate wavefunc- 
tion f(x,t) 



I. INTRODUCTION 



Two central tools in the description of a Bose- 
condensed gas are the Bogoliubov method and the 
Gross-Pitaevskii equation [0 in both its time indepen- 



dent and its time dependent guises. Because Bose con- 
densates H can now be be experimentally produced, 
there has been renewed interest in both of these tools 

EMB EM- 

As it is presently formulated, Bogoliubov's method 
treats the condensate operators a c-numbers. One conse- 
quence of this is that the resulting approximate Hamilto- 
nian does not conserve the total number of particles. The 
enforcement of number conservation in the mean only 
leads to a description which is essentially confined to a 
subspace with a single value of the mean number of par- 
ticles. 

Another difficulty is that, although there are deriva- 
tions of the validity of the time independent Gross- 
Pitaevskii equation as a description of the condensate 
ground state wavefunction (these appear largely as con- 
sequences of the adaptation of the Bogoliubov method to 
a trapped condensate), there does not appear to be any 
derivation of the validity of the time-dependent Gross- 
Pitaevskii equation as a description of the motion of a 
trapped condensate. 

This paper will solve both of these problems and show 
that their solutions are strongly connected. 

We will firstly show how to modify Bogoliubov's argu- 
ment in such a way as not to break the conservation of 
particle number. This modification yields a description 
of the particle field operator ip(x) in the form 
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V(x) « A £(x) + -=*(*) 



N 



(1) 



where A is an annihilation operator such that the eigen- 
value of A' A is N, the total number of particles, £(x) 
is the condensate wavefunction, and x(x) is a phonon 
field operator. The Bogoliubov Hamiltonian is ex- 
pressed in terms of these phonon operators, and the non- 
conservation of phonons which arises is not unexpected. 
The eigenstates of this Bogoliubov Hamiltonian are nev- 
ertheless all states with exactly N particles, so there is no 
doubt that particle number conservation is not violated. 
The method is approximate, of course, but the accuracy 
of the approximations made is exactly the same as that 
of the usual Bogoliubov method. 

It is clear that the separation of the phonon concept 
from the particle concept that enables the method to suc- 
ceed. To emphasize the validity of the the phonon con- 
cept we show that our expression for tp(x) can be used to 
demonstrate that the quantized phase of ^»(x) is the ve- 
locity potential operator, and that the operator of density 
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fluctuations is its canonical conjugate, as is well known 

The treatment of the spatially inhomogeneous case re- 
quires more care in defining exactly what the expansion 
procedure is. The Bogoliubov method and the Gross- 
Pitaevskii equations are both used to describe the weakly 
interacting Bose gas; the requirement that the gas is 
weakly interacting is formalized in our treatment by re- 
quiring the interaction potential u to be of order of mag- 
nitude 1/iV; that is, we write u = u/N, and carry out 
the asymptotic expansion in decreasing powers of ^/~N at 
fixed u. 

The treatment is then similar for the case of the time 
independent or the time dependent condensates. The 
formula (Q) is substituted in the Hamiltonian, and the 
terms of different degree in y/~N grouped together yielding 



H = NH 1 + "/NH 2 +'H 3 



(2) 



Tli is found to be a c-number, Ti.2 is linear in phonon 
operators and TI3 is quadratic in phonon operators. 

In the time-independent case the ground state is found 
by minimizing 7i 1; and this means that £(x, i) must sat- 
isfy the time-independent Gross-Pitaevskii equation — 
this also means that the linear term H.2 vanishes, and 
the excitation spectrum is given by diagonalizing the 
quadratic part H3. 

The description of the inhomogeneous condensate that 
results is almost isomorphic to that originally used by 
Fetter Q when restricted to a subspace of definite N, 
and thus is in agreement with other recent calculations. 
However, we obtain in addition a description of the rela- 
tionship between the states for different N, and in partic- 
ular are able to present the matrix elements of the field 
operators between states of N and N + 1 particles. (It 
was the need to have these in a description of condensate 
growth which motivated this work originally.) 

In the time-dependent case the logic is slightly differ- 
ent. In general there is an explicit time-dependence pro- 
portional to y/~N in x(x, t) as a result of the representa- 
tion of the time- independent V'W operator in terms of a 
time dependent condensate wavefunction £(x, t), and this 
being of order of magnitude \~N contradicts the asymp- 
totic representation in the form (0). The part propor- 
tional to \/N however can be made to cancel with the im- 
plicit time dependence arising for y/NTl^ provided £(x, t) 
satisfies the time dependent Gross-Pitaevskii equation: 
thus it is only possible to get an asymptotic expansion of 
the form ^) if £ (x, t) is a solution of the time dependent 
Gross-Pitaevskii equation. The remaining time depen- 
dence can adjusted by a suitable choice phonon mode 
functions (t) to yield a representation in terms of time 
independent phonon annihilation operators, and by going 
into an appropriate interaction picture — the condensate 
picture — the time dependence of the state vector is given 
by a Schrodinger equation H 3 (i)\$>,t) — ih~d\$,t) /dt. 
Since the number of atoms is fixed, this gives a full de- 
scription of the system; the condensate wavefunction sat- 



isfies the time-dependent Gross-Pitaevskii equation, and 
the motion of the residual quantized phonons is given by 

n 3 (t) 



II. A NUMBER CONSERVING BOGOLIUBOV 
TRANSFORMATION FOR THE 
HOMOGENEOUS BOSE GAS 

Let us start with the problem of the weakly interacting 
Bose gas confined in a large box, with no explicit trap- 
ping potential. The Hamiltonian in momentum space is 
written 



H 



k 



ki,k 2 ,k 3 ,k 4 



I ki a k 2 ak 3 a k4^k 1 +k 2 ,k 3 +k4- (3) 



The basis states are usually written 



(4) 



where n = {n^} represents the vector of occupation num- 
bers of all particles with non-zero momentum, and the 
dependence on n$ is explicitly separated. Now let us 
write them in a form which eliminates reference to n$ by 
using the total number N = no + Y]^ rik, in the form 

\N,n) (5) 

We define the operators A, Af and «k as follows 



A\N,n) = VN\N 
Af = A^A 
1 



l,n) 



Qk = 



aja k 



(6) 
(7) 

(8) 



The operators «k are essentially phonon operators, while 
A and Af are operators which refer to the total numbers of 
particles. Since the phonon operators commute with 
the total number operator Af, we can make approxima- 
tions involving the phonon operators without violating 
the conservation of total numbers of particles. 



1. Relationship between ao and A 

The operator A reduces the total number of particles 
by 1, without changing the number of k 7^ particles; 
it is thus proportional to the operator a . In fact the 
precise relationship is 



ao|JV,n) =^o|-^-l,n) 



(9) 
(10) 
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If N is very large compared with X^k^o nk — tha,t * s ' the 
system is Bose condensed — it is clear that we may make 
the approximation 



ao|^n) = A(l-^*'| \N,n) 
w A\N, n). 



(11) 
(12) 



In what follows we will be expanding in inverse powers 
of TV 1 / 2 , and it will sometimes be necessary to keep the 
more accurate form (0). 



2. Relationship between a k and Ok 
One can similarly show that one can write 



o k |JV, n) 



v^ + i-Ek^k 



N, n) 



(15) 



however we will not find it necessary to use anything 
other than the simple form (|l5|). 



3. Commutators of the Ok 

When acting on such condensed states with large N, 
we also find that we can approximate the commutator 



[ak,a£,] = <W - -^-a k a k , 



(16) 
(17) 



4- Validity of the approximations 

The approximations made are accurate to the same 
order in N as those usually assumed in the Bogoliubov 
theory. For clarity we have written equations such as 
(||,| 13) etc., as explicitly acting on a state with N par- 
ticles. However in the remainder of the paper we will 
simply write equalities between operators, which are un- 
derstood to be valid when these operators act on a vector 
in the subspace of A" particles. 



A. Transformation of the Hamiltonian 

The Hamiltonian (^) is now approximated firstly, as is 
usual in Bogoliubov theory by dropping the terms in the 
interaction which do not involve at least two operators 



ao, aj, and then secondly by making the replacements 
(and appropriate Hermitian conjugates) 



— ► TV 



ajk a k — * a k a k 



(18) 
(19) 
(20) 



in the Hamiltonian. Notice that for all of these replace- 
ments the left hand sides have essentially the same action 
on a number state as the right hand sides, but with a 
modified coefficient; that is for all of them the non-zero 
matrix elements occur for the same states, and the coef- 
ficients are almost equal, in contrast to the usual Bogoli- 
ubov transformation, which changes the essential nature 
of the operators in the number state basis. 

As in the usual Bogoliubov method, the term involving 
ajajaoao has to be treated more accuratelv^effectively by 



(13) using the more accurate approximation (|llj), giving 



alala a a = ajaoajao - a\a 



: N 2 - A/ - 2J\fJ2 a W 
k 

■ AT 2 -N- 2A/"^a k a k . 



(21) 



Although in this case it is usual to set uj — 0, for consis- 
tency with the inhomogencous case we do not make this 
assumption, and this requires a similar correction to the 
term 7ia;oao a o- The Hamiltonian ^ then becomes 

H = -u (JV 2 -JV)+Y^ - ^oKak 



k#0 



+^I]{2a k a k + a k a_ k + 4a t _ k } (22) 

k#0 

It is to be understood that this operator form is valid on 
all states, including superpositions of states with different 
numbers of particles, as long as only states with large 
eigenvalues of Af are included. The Hamiltonian can be 
diagonalized for any eigenvalue of A/", to get the usual 
Bogoliubov spectrum; the eigenstates are simultaneous 
eigenstates of Af and H. 



1. Asymptotic expansion 

The expression of the Hamiltonian in the form ( p2f ) can 
be put in the form of an asymptotic expansion in inverse 
powers of A" in the sense that the restriction of H to a 
subspace of definite A~ is approximated by (^||) for suffi- 
ciently large N, provided it is understood that u is itself 
a small quantity of order of magnitude 1/N . Under that 
condition it can be checked that approximations made 
in its derivation affect only the terms of lower order in 
A^; i.e., the corrections are of order 1/N. The require- 
ment that u be "small" amounts to an assumption that 
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the kinetic plus trapping potential terms, and interaction 
energy terms in (|22|), are of the same order of magnitude. 
A formal statement of this requirement can be made by 
writing u = u/N in (^2|), and then solving the problem by 
making an expansion in appropriate decreasing powers of 
N. 

The Hamiltonian in this form has two major advan- 
tages over the deceptively similar conventional form. The 
most significant is that the conservation of total numbers 
of particles is maintained; the quasiparticles on the other 
hand appear only as phonons, and are not conserved. 
The second major advantage is that the approximation 
method is a systematic expansion in inverse powers of 
the large quantity N, which has a definite value in any 
case one considers. 



2. Expression of field operators in terms of the velocity 
potential operator 

The consequences of this method on the representation 
of the field operators are interesting. We can show that 
the representation of phase and density fluctuations as in 
H is a natural consequence. 

The field operator is 



^00 = 77^ I "n 



A_ 
V 



-y 

k/0 



a k e 



ik-x 



(23) 



(24) 



We now introduce the Bogoliubov transformation, ap- 
proximated by assuming k is very small (though this is 
not an essential assumption, and is introduced only to 
make it match up with ||), as 



hk 



y/2hk/mv(N) y 2mv(N) 



(25) 



where v(N) = y/ uN / m is the speed of sound for long 
wavelengths. Clearly, the second term is much smaller 
than the first for small k, so we keep only the first part 
initially. We then find that 



V-(x) « -= 



A_ 

V 



k^O 



mv 
2hkN 



& k e" lkx ) 



-= exp (i$(x)) 
v V 



(26) 
(27) 



since 



k^O 



is small. 

This matches up exactly with (27. f ) of ||. That is, the 
operator $ is the quantized velocity potential operator, 
with the substitution 6 k — * ic^ to match up with their 
notation. If we now include the next order term we find 
density fluctuations as well. Using the full form (p5|) for 
the density 



N 



V V 



Po + Sp(x) 



1 T [NM r 
V ^ V 2mv I 

k#0 



kx 



t c -ik-x I 



(29) 
(30) 



giving the correct operator for the density fluctuations, 
as in d (24.10). 

In this limit that only long wavelengths are involved, 
the Hamiltonian ( p2| ) can be rewritten in terms of a part 
related to sound waves, as in ||, and an additional purely 
J\f dependent part: 



1 

H = -u 
2 



2 ' 



d 3 x 



p (V.$(x)) 2 + v (AA) 2 ^ 

Po 



(31) 



Here Eq (Af) is the Bogoliubov ground state energy. 

Notice that not only do po and v(N) depend on N, but 
that the very definitions of the $(x) and 5p(x) also de- 
pend on N. However, if we believe that a superselection 
rule applies whereby only eigenstates of N occur in na- 
ture, the correct treatment of states with uncertain total 
numbers of particles must arise from an incoherent su- 
perposition of solutions, for various N, of the equations 
of motion arising from (|3l|). 

III. THE SPATIALLY INHOMOGENEOUS CASE 

The spatially inhomogeneous situation, which can arise 
either by the existence of vortices or because there is a 
trapping potential confining the gas, was first treated by 
Fetter 0, and in the case of a trapped gas, has been 
treated more recently in j|,[Il]-[l3|] . 

The standard formalism, as used by Fetter JtJ , is based 
on the work of Hugenholtz and Pines p6| which in sum- 
mary, shows that the correct ground state and excitation 
spectrum is given by the following prescription: 



1. Replace a 



N, a 



N in both H and Af. 



2. For a given p, find the value of N and the state 
which minimize (K), where K = H — pN . This 
gives the ground state energy. 

3. The excited states are given by the higher eigen- 
states of K with aj — > VN, ao — > y/~N. By neglect- 
ing terms in N of order less than N, one obtains 
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the Bogoliubov Hamiltonian, which can be exactly 
diagonalized, but Hugenholtz and Pines in fact in- 
cluded more terms than these in their evaluation of 
the ground state energy, and hence obtain a more 
accurate result. 

It is particularly important to note that there is no basis 
for using K without setting — > y/~N, ao — > y/N, as was 
done by [p"l| , which can be seen as an initial attempt to 
get a Bogoliubov description which covers a range of N 
values. 



A. Formulation of the modified Bogoliubov method 

1. Expression of the field operators in terms of phonon 
operators 

Our treatment will be based on that already used for 
the spatially homogeneous case. We therefore consider 
the general case for which the Hamiltonian is 

H=-— [ d 3 xV> f (x)V 2 VXx) + / d 3 x^(x)l/(x)(/)(x) 
2m J J 

(32) 

We will make the substitution 



+- / d 3 xV' t (x)V' t (x)V'(x)V'(x) 



^(x) = a 




(33) 
(34) 
(35) 



Here the operators {ao, at} are independent creation and 
destruction operators satisfying the usual creation and 
destruction operator commutation relations, and the set 
of functions {£o( x ), £fc( x )} is a complete orthonormal set. 
This is necessary to produce the mandatory field operator 
commutation relation [^(x), V^(x')] = S(x — x'). 

The meaning of A and ctk is essentially as previously 
defined except that the exact nature of the states involved 
is not yet defined. The basic issue is that the occupation 
of the mode described by ao is very large compared to 
the occupation of any of the modes described by the a*;. 
Under these conditions, apart from the notation change 
k (momentum) — > k (an arbitrary label) the relationship 
between the operators A, ao, a/., a-k is exactly the same as 
in Sect. 

Notice that X (x) has a non-local commutation relation 

[ X (x), X t(x')]=]>>(x)a(x') (36 ) 

k 

= <5(x-x')-£(x)r(x') (37 ) 

= i?(x,x'), (38) 



which arises because the X operators act only in the 
subspace orthogonal to the condensate wavefunction £. 
These are of the same form as those for similar operators 
considered by Fetter M . 



2. Expansion of the Hamiltonian 

The Bogoliubov method will be valid for large N and 
small u, and, as noted above, this is expressed more pre- 
cisely by setting u = u/N, and expanding the Hamil- 
tonian in decreasing powers of N for fixed u, after sub- 
stituting for the field operators using (|34|). In order to 
eliminate ajao as in the homogeneous case, we note that 



J2 a k a k = J rf 3 xx f (x)x(x) 



(39) 



Carrying out this procedure, we then get 

H = NH1 + VNH2+H3 + ... (40) 

in which 

U x = ~ J d 3 xr(x)V 2 £(x) + J d 3 xr(x)F(x)C(x) 

+ ^ J d 3 x|£(x)| 4 , (41) 
n 2 = ~ [ d 3 x{ X t(x)V 2 £(x) +r(x)V 2 X (x)} 

+ J d 3 xy(x){ x t( x ) ? ( x ) +r(x ) x ( x) } 

+u J d 3 x{|e(x)| 2 ^(x)x t (x) + |e(x)| 2 r(x)x(x)}, 



(42) 



Hz = - 



2m 



J d 3 x X t( x)V 2 x (x)+ J d 3 x X t (x)U(x) X (x) 



+ yd 3 x|^(e(x) X t(x)) 2 + ^(r(x)x(x)) 2 
+ X + (x)x(x)(2fi|S(x)| a - A i)|-| f d 3 y|£(y) 



(43) 



in which 



n 2 v 2 

2m 



C(y) 



d 3 y(c(y)( 
+v(y)\ay)\ 2 + u\C(y)\ 4 ) 



(44) 



If we now choose ^ to minimize 7i±, subject to the con- 
dition that / d 3 x |£(x)| 2 = 1, we find that 

1. The condition for a local minimum of 7ii is that 
the condensate wavefunction £ (x) satisfies the time 
independent Gross- Pitaevskii equation 



5 



2m 



V 2 £(x) + F(x)£(x) + fi|£(x)|^(x) = mC(x). 

(45) 



Here \x arises as a Lagrange multiplier necessary 
to maintain the normalization of £(x) — any valid 
value of fi is permitted in this procedure, and the 
value obtained is consistent with (44) 



2. Under this condition the terms linear in \ or 
vanish, since the £fe(x) are a se t of functions or- 
thogonal to £(x). 

Since TV is known, p5h is to be considered as a nonlinear 
eigenvalue equation for £(x). Thus, the possible values 
of // are the eigenvalues of the Gross-Pitaevskii equation. 

For stability 7i 3 should be positive definite; whether 
this can be satisfied seems to depend on the actual so- 
lution of the Gross-Pitaevskii equation which is under 
consideration. If the condition is not satisfied, the state 
represented by the solution of the Gross-Pitaevskii equa- 
tion is not stable, and the extremum found is not a min- 
imum. 

When restricted to a fixed TV subspace, the expansion 

f,pl|-^3|), is exactly the same as that obtained by Fetter 
using the method based on the minimization of (K). 



3. Diagonalization ofTiz 

The expression (^) for TL 3 can be diagonalized most 
simply by working in the basis, in much the same way 
as Javanainen This method is also preferred here 

since it generalizes rather straightforwardly to the time- 
dependent situation. Thus one can write 



n 3 = E 3 



k,q I 



kq + Fkq) a\a q 



where 

Lkq 



d 3 xS(x)^-v 2 + y(x)K,(x) 



F kq =uj d 3 x£*(x) {2|C(x)| 2 - fi} e,(x) 
u J d 3 x^(x){£(x) 2 }£ 9 (x) 



Gkq 
E3 



d 3 y\ay)\ 4 



(46) 

(47) 
(48) 
(49) 
(50) 



4- Quasiparticles 

The Hamiltonian H3 can be diagonalized by a Bogoli- 
ubov transformation of the form 



Oik = E C kmb m + E S kmb\ n 



(51) 



and here b m is a quasiparticle destruction operator. We 
can then write 

= (P™(*)b m + q m {*)bl) (52) 

m 

with 

PmW = ^ c fc™&( x ) ( 53 ) 
k 

q m (x) = ^ Skm^k (x) . (54) 
k 

The diagonalized Hamiltonian is then written 

H 3 = hu g (N) + he m (N)blb m . (55) 



Notice that, even though the TV dependence is not al- 
ways explicitly written, almost everything in the above 
is a function of TV. 



B. Relationship between ground states for TV and 
N + 1 particles 

The operator A contains an TV dependence which arises 
from the change in the shape of the ground state wave- 
function £(x) as TV changes. This is an effect which does 
not arise in a spatially homogeneous situation. Thus we 
can exhibit this feature by writing, as a result of (p3), 



a (N) 



d 3 xr(x,Aov(x), 



(56) 



where the explicit dependence on TV of both £ and clq has 
now been written. 

This means that aJ(TV)|TV, n = 0) is a state with TV+ 1 
particles in the wavefunction corresponding to ground 
state of the TV particle state; it is not the ground state 
for TV + 1 particles. In other words, aJ(TV) is not the 
operator which converts an TV particle ground state into 
an TV + 1 particle ground state. 

We will therefore compute the appropriate operator. 
We can write 

ao(JV + l) =a (N) 

+ yd 3 x(r(x,TV + l)-r(x,iV))^(x) (57) 

j3 ar(x,Tv) 



ao(N) 



d J x- 



<9TV 



-V>(x). 



(58) 
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Now expand <9£* (x, N) /dN in the N particle basis states 
as 



ag*(x, N) 
dN 



k 



r k e k (x,N)} (59) 



The requirement that £ be normalized allows us to 
deduce that ro is real, and by redefining £(N) — > 
£(iV) exp{— i JdNro(N)/N} we can transform ro to 
zero; that is, we can choose a family of solutions of the 
time-independent Gross-Pitaevskii equation in which the 
relative phases of the members of the family are such 
that ro = 0. Once ro has been eliminated, it is clear that 
r k are of order of magnitude 1. This means that we can 
write approximately 



a (N + 1) « a (N) + i ^ r fc a fc (JV) 



(60) 



We now want to relate the N particle ground state to the 
N + 1 particle ground state. We can write 



| AT) 



JV 



-|0) 



|JV + 1>jv+i = 



y/Nl 
{al(N + l)} 



N+l 



(61) 



(62) 



We first write 
a k (N) thus 



(J) in terms of the phonon operators 



a (N + l) =a (N) 1 + 



£ 



r k a k {N) 
N Jn 



(63) 



and raise this to the power N + l, using the binomial 
theorem to get, accurate to order 1/N , 



{<4(AT+l)} 



N+l 



k 4(N) 



N 



N+l 



(64) 



From (pUB2n this means that 



r i N+1 r i " 

\al(N + l)\ U(N)\ 

\N + 1) N+1 = ± L 1 / — ■ |0) (65) 



'N\ 



1 + 



E 



4w 

VnTi 



\N) 



N 



(66) 



Similar procedures can be carried out for the a k (N) and 
hence the a k (N) , but the differences between these at N 
and N + l will have a negligible effect because of their 
small occupation numbers. 

From the above, we can say that the operator B'(N) 
which connects the N and N+l ground states through 



B^(N)\N) N = VN+~T\N + 1) N+1 (67) 
is given approximately by (to order 1/N) 



B*(N) 



E^4 )4m 

k ) 



(68) 



This means that the field operator expansion ( p5| ) now 
takes the form 



-•(xl ^ B{X) (e(x) + -^Xs(x) 



where 



Xb(x) = x(x) -C(x)^ 



and 



= ^2 (/™( X ) 6 ™ + .9m(*)&L 

m 

/ m (x) = J?m(x) - $(x) ^ r ^kra 
m 

5m(x) = <? m (x) - C(x) r fc s fem 



(69) 

(70) 
(71) 

(72) 
(73) 



The physical distinction between the quasiparticle wave- 
functions p m , q m and the modified wavefunctions f m , gm 
is that the first set are amplitudes for the production or 
absorption of quasiparticle without changing N, e.g., by 
application of sound waves, whereas the second set give 
the corresponding amplitudes for changing quasiparticle 
numbers by changing the particle number by 1. 

It is particularly interesting to see that the corrections 
to the quasiparticle term are of the same order of mag- 
nitude as the original terms; the correction is thus very 
significant. It reflects the fact that the change in the 
groundstate wavefunction from N to N + 1 particles af- 
fects the N particles already present as well as the added 
particle. 



IV. TIME-DEPENDENT SOLUTIONS 



Suppose we now consider a representation like (|33J-|35j) 
in which however we have a time-dependent condensate 
wavefunction £ (x, t) . This is with fixed N since we are 
not considering condensate growth; rather, we are in- 
vestigating the situation where the condensate has been 
macroscopically disturbed from the stationary state as in 
the recent experiments in JILA [Q and MIT Q. 

Let us proceed within the Schrodinger picture as fol- 
lows. Firstly suppose the condensate wavefunction (nor- 
malized to 1) is £(x, t) and define condensate and non- 
condensate parts of the field operator by 
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lfo(x,t)=£(x,t) / d 3 x'r(x',t)V(x') 



= a (^(x, «) 

Ipncfct) = 1p(x) - ^ (x,t). 



(74) 

(75) 
(76) 



We can set up non-condensate mode operators by defin- 



ing 



a k (t) = / d 3 x^(x,t)V(x) 



(77) 



where {£fe(x, £)} are any orthonormal set of wavefunc- 
tions orthogonal to the condensate wavefunction £(x, t). 
The total number operator is 



TV = / (ftc^x^x) 



(78) 



and is time independent. We can correspondingly define 
condensate and non-condensate number operators by 



d 3 xtl)l c (x,t)tp nc (x, t). 



(79) 
(80) 

(81) 



The time-dependence of these operators arises because 
their definitions, which involve the time-dependent wave- 
functions £(x, t), £fc(x, i), change with time; they are how- 
ever still operators in the Schrodinger picture. 

Using these definitions it is possible to define the op- 
erator A essentially as before so that, 



a (t) 



N DC (t) 
N 



A(t). 



(82) 



Although A(t) is time-dependent, N = A'(t)A(t) is time 
independent. The time-dependence of A{t) comes about 
because of the changing definition of the modes unaf- 
fected by its action. That is, A reduces the total number 
of particles while leaving the number of non-condensate 
particles the same, but the definitions of the condensate 
and non-condensate modes are themselves time depen- 
dent. Further, although the operators for condensate and 
non-condensate modes commute, the A(t) operator does 
not commute exactly with any of ao(i), %(£), (Zjfe(i) or 

4(t). 

It will be particularly important to keep track of all 
time-dependences which are a result of the the expres- 
sion of the time-independent Schrodinger picture field 
operators ^>(x) as projections on the time-dependent ba- 
sis vectors {£(x, i), £fe(x, i)} — the only operator in this 
basis that does not develop such a time-dependence is 
the operator Af. We will show that a unitary transfor- 
mation can be introduced which transforms us to a kind 
of interaction picture, which we shall call the conden- 
sate picture, in which we have a description in terms of 
time-independent phonon creation and destruction oper- 
ators /3fc,/3fc corresponding to the time-dependent modes 



£/c(x, t) provided that the condensate wavefunction £(x, t) 
satisfies the time- dependent Gross- Pitaevskii equation 

^§ = -^V 2 £ + n + €|£| 2 £ + (83) 

where p,(t) is arbitrary, and of course can be elimi- 
nated by multiplying £ by an appropriate time-dependent 
phase. 



A. Time dependence of the phonon operators 

1. Explicit time-dependence 

Since {£(x, i), £fe(x, t)} are a complete orthonormal set, 
we can define a projector onto the non-condensate modes 
by 



R(x,x',t)=j2^,m(^',t) 

k 

= J(x-x')-£(x,i)£V,*) 
so that, using a time-dependent version of (pi 



A(i)-^=x(M)« / d 3 x'R(x,x',t)i>(x') 



(84) 
(85) 

(86) 



and this means that to the lowest order in v N we can 
write 



<(x, *) = --!= ^ (t) J d 3 x'i*(x,x',i)V(x') 



(87) 



The explicit time-dependence of \ arises from the pro- 
jector and the operator A'(t). We use expression 



i4t(t)woS(t)= / d 3 x^(x,t)^(x) 



so that 



8A\t) 
dt 



d 3 x^(x,t)V' t (x) 
At(t) / d 3 xe(x,<)r(x,<). 



(88) 

(89) 
(90) 



Now using the explicit form (|35|) for the projector R, and 
resubstituting for ij> in terms of A and x we deduce 



dVe(x',t)r(x',i) x(x,*) 



dt 

+ J d 3 x'7?(x,x',<){\/iVC(x',t) + x(x',<)} 

(91) 

This time-dependence is only the explicit time- 
dependence, and does not include that arising from the 
commutator with the Hamiltonian 
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2. Expansion of the Hamiltonian and the condensate picture 

The Hamiltonian can now be expansied in decreas- 
ing powers of y/~N in exactly the same way as for the 
time independent case, with appropriate substitutions. 
However, the largest terms arising from ( |42| ) and j9l| ) 
are both proportional to ^/N, and these must be ar- 
ranged to cancel, since otherwise x/VN would develop a 
c- number component comparable with £, invalidating the 
expansion ( [35]). We will show that this cancellation is 
only possible if £ (x, i) satisfies the time-dependent Gross- 
Pitaevskii equation. 



The coefficient of y/N in (91) can be written as 



j d 3 x'i?(x,x',t)£(x',i) = % -[H R (tlx{*M (92) 



in which 



H R (t)=iH ( d 3 x{e(x,<) X t( Xj< )_^( x ,t) x ( Xji )| 



Here we have used the identity 



(93) 



= J d 3 x'i?(x,x',i)£(x',i) + / d 3 x'i?(x,x',i)£(x',i), 

(94) 

which arises from the fact that / d 3 x'i?(x, x', t)£(x', t) = 

o.) 

If £ satisfies the time-dependent Gross-Pitaevskii equa- 
tion (^3|) , it is clear that TLr = —H.2 ■ This means that we 
can make a unitary transformation from the Schrodinger 
picture to a picture which we call the condensate picture, 
defined by 



\$,t)^V(t)\$,t) = \$,t) c 
^(x)^F(i)V(x)F- 1 (t) = ^ c (x,t) 



dt 



V(t) 



NH R (t)V(t). 



(95) 
(96) 

(97) 



In this picture the part of the explicit time-dependence 
of x proportional to \/~N in (|9l]) will disappear, and the 
term \/~N H2 in the equation of motion for the states 
|$,t) c will cancel with a term arising from the unitary 
transformation. We will then be able to write the equa- 
tion of motion in this picture as (to order N ) 



(NHi+Hs(t))\^,t) c = ih-\^,t) c 



in which 



w i = "^ / d 3 xr(x,t)V 2 £(x,t) 



+ / d 3 xnx,^(x)£(x,i) + ^ / d 3 x|£(x,i)| 4 , (99) 



H 3 (t) = Jd 3 

+ X t(x,i)F(x) Xc (x, t) 



— X t(x,i)V 2 Xc (x,i) 



+ |(e(x,i)^(x,t)) 2 + ^(r(x,t) Xc (x,t)) 2 



+xl (x, <) Xc (x, t) [2fi|£(x, t)\ 2 - IM(tj\ 



u 2 !*y\tM\* 



(100) 



in which fj,(t) has no connection with the /i(i) in 
but is given by 



K*) = j rf 3 y(r(y,t)(- 

-'I' I d 3 yC(y,t) 



2m 



£(y,t) 



(101) 
(102) 



Notice that no time-dependence is written for Hi, since 
it is in fact a constant c-number when £ satisfies the time- 
dependent Gross-Pitaevskii equation. 

Note that although the equation of motion is given 
by (|98|), the energy is given by the full Hamiltonian 
H = NTL\ + v^/V H2 + H%(t), with the operators evalu- 
ated in the condensate picture. Since there is no time- 
dependence of the Schrodinger picture Hamiltonian, the 
total energy must be conserved. 



3. Time-independent phonon operators 

The expansion of the phonon field in the condensate 
picture as 



X. 



: (x,i) =£)&(x,t) 



OLk 



(103) 



does not automatically require that ak be time indepen- 
dent. However the mode functions have been so far 
essentially arbitrary. The explicit time-dependence of Xc 
in the condensate picture is given by omitting the term 
prop ortion al to \/N in (|^), and by inserting the expan- 
sion ( |103[ ) into this, we find that can be chosen to 
be independent of time if the mode functions satisfy the 
equation of motion 



a&(x,t) 
dt 



£(x,t) / d 3 yt(y,t)Uy,t) 



(104) 



It is straightforward to check that the requirements that 
£fc form an orthonormal set, and are orthogonal to £ are 
both preserved by this equation of motion. 
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S ubst ituting the expansion ( |103| ) into Ti 3 {t) as given 
by (100) we obtain 7is(t) in the form 



H 3 (t) = E 3 (t) + ]T { (L kQ + F kq (t)) 4« 

k,q 



+G kq {t)a\a\ + G* kq (t)a k a 



(105) 



Here Lk q (t), Fk q {tV Gk q (t) and E${t) are defined in the 
same way as in (pf7|-p0[), but using in this case the con- 
densate picture operators and the time-dependent mode 
functions which are now arbitr ary o nly at the initial time, 
since the equation of motion (104) gives them for all fu- 
ture times. 



B. Summary of the time-dependent description 

We have shown that the field operators can be writ- 
ten in a time-dependent version of the form (|35"|). This 
expansion is valid as an asymptotic expansion in 
if and only if £ (x, t) satisfied the time dependent Gross- 
Pitaevskii equation. 



V. CONCLUSION 

The adapted Bogoliubov method presented in this pa- 
per gives a precise meaning to the definition of the macro- 
scopic wavefunction usually defined by the limiting pro- 
cedure fl3l 



(N,m,t\ipi(x)\N- l,m,t) ~ y/N£*(x,t), 



(106) 



where, in the words of ]3|], 

"...since the condensate contains a macro- 
scopically large number of particles, changing 
this number by 1 does not essentially affect 
the state of the system; we may say that the 
result of adding (or removing) one particle in 
the condensate is to convert a state of the 
system of N particles into the 'same' state of 
a system of N + 1 particles. ... the symbols 
\N,m,t) and \N + l,m,t) denote two 'like' 
states differing only as regards the number of 
particles in the system." 

This definition has the disadvantage of being rather 
vague about exactly what is meant by the concept of the 
'same' or 'like' states which differ only in that the values 
of N differ by 1. From this paper it is clear that the 
description of the eigenfunctions of the non-condensed 
particles depends on the value of N, and although the 
difference between the eigenfunctions for N and N + 1 
must become negligible for N —* oo, if we are contem- 
plating a situation in which the condensate is growing, 



as in recent experiments, this difference may have non- 
negligible effects. 

Our method can also be seen as the logical completion 
of the aim of Hugenholtz and Pines j2(| to treat the Bose 
condensed gas by elimination the ground state. Their 
method was only able to do this at the expense of break- 
ing the exact conservation of particle numbers, whereas 
our method puts the approximation in the relationship 
between the operators A, a k and the particle operators, 
and in their exact commutation relations. However, there 
is no good reason to why one should not in principle 
seek a more accurate approximation than has been done 
here by computing terms of higher order in the inverse 
VN expansion, and get results comparable with those of 
Hugenholtz and Pines. For the condensates at present in 
existence this is not an urgent problem; rather the main 
problem is to get as simple a description as possible of 
the eigenfunctions as a function of N as well as the other 
variables. 

Griffin [Q has recently shown how the Hartree-Fock- 
Bogoliubov method can also be used. This is closely 
related to any Bogoliubov method, including this one, 
but it suffers from not being a systematic expansion in 
any parameter. For example, depending on assumptions 
made on certain averages one may or may not obtain a 
gapless spectrum. A systematic method would identify 
all terms of a given order in an appropriate small pa- 
rameter, and either use all of them or use none of them. 
If the Hartree-Fock-Bogoliubov method is developed in 
inverse powers of VN, one will find the same results as 
presented here to the degree of accuracy presented here. 

It is well known that the Bogoliubov method is gapless 
p7| , p8[ (in the long wavelength limit, the energy levels 
approach the energy of the ground state; i.e., a phonon 
of very long wavelength has vanishingly small energy). 
The method of derivation used here shows that the non- 
conservation of particle numbers, often seen as the hall- 
mark of the method, is not at all essential to the method; 
that it is really takes only a slightly different point of 
view to see that the non conservation of particle num- 
bers arises from a slightly inappropriate way of looking 
at the quasiparticle operators. Nevertheless, the elimi- 
nation of this problem is absolutely essential if we wish 
to study condensate growth. Approximations which do 
not preserve particle number conservation can give rise 
to spurious terms in the equations of motion for the de- 
velopment of a condensate, and thus make it difficult to 
identify the true details of the growth process. 

The reason for developing the modified Bogoliubov 
method is thus to apply the method to the growth of 
the condensate, and this will require its incorporation 
into the framework of Quantum Kinetic Theory ]3lf| . If 
the condensate growth is rather slow, it will adequate to 
use the time-independent formalism of Sect. [II, but for 
faster condensate growth, the time-dependent formalism 
will be indispensable. These aspects will be treated in 
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